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We briefly discuss the quark-antiquark Bethe-Salpeter equation and the 
quark Dyson-Schwinger equation derived in preceding papers. We also con- 
sider the qq quadratic mass operator M 2 = {w\ + W2) 2 + U obtained by three- 
dimensional reduction of the BS equation and the related approximate center 
of mass Hamiltonian or linear mass operator Hqm = M = W1+W2 + V + ■ ■ ■■ 

\Q \ We revue previous results on the spectrum and the Regge trajectories ob- 

tained by an approximate diagonalization of Hqm an d report new results 

{Sj ' similarly obtained for the original M 2 . We show that in both cases we suc- 

ceed to reproduce fairly well the entire meson spectrum in the cases in which 

q^ . the numerical calculations were actually practicable and with the exception 

of the light pseudoscalar states (related to the chiral symmetry problematic) . 
A small rearrangement of the parameters and the use of a running coupling 
Q_ ( ' constant is necessary in the M 2 case. 
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I. INTRODUCTION 

In preceding papers |T|Q] (cf. also |3[|), using a Feynmann-Schwinger type of path integral 
representation and an appropriate ansatz for the Wilson loop correlator, we have obtained a 
Bethe-Salpeter and a Dyson- Schwinger equation for certain "second order" quark- ant iquark 
and single quark Green functions. 

To solve the qq bound state problem one should in principle solve the DS equation and 
use the resulting quark propagator in the BS equation. In practice even an approximate 
treatment of such problem in its full four-dimensional form seems to be extremely hard and 
for we have to resort to the use of free propagator and of a three-dimensional reduction 
of the BS equation (instantaneous approximation). Such reduction takes the form of the 
eigenvalue equation for an effective squared mass operator M 2 = {w\ + W2) 2 + U, Wi and 
W2 being the relativistic free energies of the quarks and U an interaction related to the BS 
kernel. 

In more conventional terms one can also consider a center of mass Hamiltonian Hqm = 
M = W\ + W2 + V + . . . , where at the lowest order V differs from U only for kinematic 
factors. This last form can be more directly compared with usual relativistic and non 
relativistic potential models and V turns out to have various significant limit expressions. 
In the static limit V takes the Cornell form 

V = -\^ + ar. (1.1) 
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In the heavy masses limit, by an — expansion (and an appropriate Foldy-Wouthuysen trans- 
formation) it reproduces the semi relativistic potential discussed in || and @ . If the spin 
dependent terms are neglected, it becomes identical (apart from a question of ordering) to 
the potential corresponding to the relativistic flux tube model 0, up to the first order in 
the coupling constant a s and the string tension a. 

In || we have solved numerically the eigenvalue equation for Hqm, as we shall explain 
later, neglecting the spin-orbit terms but including the hyperfine separation. We have suc- 
ceeded to reproduce fairly well the entire heavy-heavy, light-light and light-heavy quarko- 
nium spectrum and Regge trajectories when the actual calculations were feasible. The only 
real exception was the case of the ground light pseudoscalar mesons, for which the three- 
dimensional reduction of the BS equation does not seem to be appropriate, due to the chiral 
symmetry breaking problem. 

Concerning the choice of the constants, the light quark masses were fixed on typical 
current values, m u = m s = 10 MeV, m s = 200 MeV and only the heavy quark masses, 
the strong coupling constant and the string tension were used as fitting parameters. Good 
agreement with the data was found for m c = 1.40 GeV, m& = 4.81 GeV, a s = 0.363, a = 
0.175GeV 2 . 

In spite of the success attained it turns out,however, that the quantity (V^ 2 ) is not 
negligible bringing e.g. to corrections ranging between few tens and 150 MeV in the cc 
case. For this reason in this paper we have repeated the calculations for the more complex 
operator M 2 . A good agreement is again obtained at the price of a small rearrangement of 
the parameters and of using a running coupling constant given by the usual perturbative 
expression 



"■«>= (11 -|N f)ln g ™ 

cut at a maximum value a s (0). 

We have taken N t = 4, A = 200 MeV, a 8 (0) = 0.35 and a = 0.2 GeV 2 , where the last 
two values have been chosen in order to reproduce the correct J/^> — r\ c separation and the 
Regge trajectory slope. We have also chosen m c = 1.394 GeV and m& = 4.763 GeV in order 
to reproduce exactly the masses of J/\l/ and T; on the contrary we have left unchanged the 
masses of the light quarks. 

In the following we revue the BS and the DS equations in sect. II and the three- 
dimensional reduction in sect. III. In sect. IV we report and discuss the results obtained in 
|| for H CM and the new results for M 2 . 



II. BETHE-SALPETER AND DYSON-SCHWINGER EQUATIONS 

The gauge invariant "second order" Green functions considered in HJ, @ and |J were 
defined as 

H*(xi,x 2 ;yi,V2) = -~Tr c {U{x 2 ,x 1 )A c [{x l ,y 1 ;A)U{y 1 ,y 2 )A°{x 2 ,y 2 ;-A)), (2.1) 

Hz\x -y) = iTr c (U(y, x)A°(x, y; A)) , (2.2) 

where the tilde and Trc denote the transposition and the trace, respectively, over the color 
indices alone; U is a path-ordered gauge string joining a to b (Schwinger string), U(b,a) = 
P exp lig J a dx^ A^{x) [; while A a (x, y; A) stands for the "second order" quark propagator 
in a external gauge field A^, defined by the iterated Dirac equation 

(£>„£>" + m 2 - l -g o^F» v )A°(x, y; A) = -5\x - y) , (2.3) 

with a^ u = |[7 M ,7^]; finally the angle brackets in Q2.1| ) and ( |2.2| ) denote average on the 
gauge variables alone (weighted in principle with the determinant Mf(A) resulting from the 
explicit integration of the fermionic fields) . 

The quantities H sl (xi,x 2 ;yi,y 2 ) and H gl (x — y) are simply related to their ordinary 
"first order" counterparts and can be equivalently used for the determination of the bound 
state. The advantage they offer is that of admitting path integral representations in terms 
of quark world lines which derive from the similar Feynmann- Schwinger representation for 
A a (x, y\ A). Such representations depend on the gauge field only trough Wilson correlators 
W[T] = ^{TrP exp{ig §r dx^Afj}) , associated to loops T made by the quark and antiquark 
world lines closed by "Schwinger strings". 

In principle, as a consequence, the above correlators should determine the whole dynam- 
ics. Unfortunately, due to confinement and the consequent failure of a purely perturbative 
approach, a consistent analytic evaluation of W from the Lagrangian alone is not possi- 
ble today and one has to rely on models based on incomplete theoretical arguments and 



lattice simulation information. The most naive, but at the same time less arbitrary assump- 
tion, consists in writing i In W as the sum of its perturbative expression and an area term 
(modified area law (MAL) model) 

t\nW = t(\nW) pcvt + aS min , (2.4) 

where the first quantity is supposed to give correctly the short range limit, the second 
the long range one. Notice in principle any more sophisticated model could be used, at 
the condition it preserves certain general properties of functional derivability of the exact 
expression. In practice not even ( |2.4| ) can be treated exactly. Actually one has to replace the 
minimal surface S m in by its "equal time straight line approximation" , defined as the surface 
spanned by straight lines joining equal time opposite points of the loop T in the qq center 
of mass frame. 

The path integral representations obtained in such a way could be used directly for 
numerical calculations or for analytic developments. In the last context it is convenient to 
consider a second type of second order functions H(xi,X2,yi,y2) and H(x — y) obtained 
from H gl (xi,X2,yi,y2) and H gl (x — y) by omitting in their path integral representation the 
contributions to i In W coming from gluon lines or straight lines involving points of the 
Schwinger strings. In the limit of vanishing X\ — X2, y\ — 2/2 or x — y such new quantities 
coincide with the original ones and are completely equivalent to them for what concerns the 
determination of bound states, effective masses, quark condensates, etc. 

By an appropriate recurrence method, an inhomogeneous Bethe-Salpeter equation and 
a Dyson-Schwinger equation can be derived for H(x\, x 2 ,yi, 2/2) and H(x — y), respectively. 
In the momentum space, the corresponding homogeneous BS-equation can be written (in a 
4x4 matrix representation) 

_ ,, , . f d^u f . , 1 k + u 1 „ k + u. 

#i(ip + k)o"$p(u)a b H 2 ( - l -P + k), (2.5) 

where $p(ft) denotes an appropriate wave function and the center of mass frame has to be 
understood; i.e. P = (m B ,0), m B being the bound state mass. Similarly, in terms of the 
irreducible self-energy defined by H{k) = H (k) +iH (k)T(k)H(k) the DS-equation can be 
written also 

f (*) = j^yM k - z - ^T' ^-ym ° b ■ (2.6) 



Notice that in principle ( |2.5| ) and (|2.6|) are exact equations. However the kernels I a b 
are generated in the form of an expansion in a s and a. At the lowest order in both such 
constants, we have explicitly 

Io;o(Q;p,p) = l67i-a s p a p >f3 D af3 (Q) + 

+4a J d*Ce- lCl <\CHpoW ) jf 1 d\{plp' 2 - [Ap' p T + (1 - AW T ]¥ 



4 
I^{Q]P,p') = 47ri-a B (<*£Q„ - S^Q^D^Q) - 

U, P a(Q]P,p') = -Am^a s p a (5^ p Q a - S^Q p )D aP (Q) + 
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Iw P cr(Q;p,p') = ^a s (6«Q v - S:Q,)(S a p Q a - 5 a a Q p )D a ^Q) (2.7) 

where in the second and in the third equation Co = has to be understood. Notice that the 
use of (|1.2Q in (|2.7|) would amount to include higher order contributions. 



III. THREE-DIMENSIONAL REDUCTION OF THE BS EQUATION 



To obtain from (|2.5| ) a three-dimensional equation we can perform on such equation the 
so called instantaneous approximation. This consists in replacing in ( |2.5|) H% (p) with the 
free quark propagator 2 ~^ 2 and the kernel I a b with J^ st (k, k') obtained from I a b setting 

ko = k' = r] 2 Wl 2 Wl — iji " 2 2 ^ 2 with iUj = Jm 2 + k 2 and iw^ = Jm 2 + k' 2 . Then, by 
performing explicitly the integration over k' and further integrating the resulting expression 
in fc , we obtain 

(u7l + W 2 )Vm S ( k ) + 



+ / THjih; 4b ( k > k h/ o / / a VW 1 *)^ = m B¥>m B ( k ), (3-1; 

' (27r) ci y 2w\W2 V 2w[w 2 



with ^p(k) = v /&/- oo ^ $ P (/c). 

Eq. ( p.l|) is the eigenvalue equation for the squared mass operator, 



M 1 = M' + U (3.2) 



with M = Jm\ + k 2 + Jm\ + k 2 and 



1 / Wl + W 2 finst ,, , a WJ+W2 _g_b 

(2tt) 3 V 2 Wl «; 2 ab l ' J V 2u« 



/virrlLA _ x 1 W2 pMt/i, ua w i T_ W2 ~"~fc cq ^ 

<W0 = 7S3s\/-o^r-^ ( k ' k )V^T^77- a i a 2- (3-3) 



The quadratic form of Eq. (|3.2|) obviously derives from the second order character of the 
formalism we have used. 

In more usual terms one can also write 

H CM = M = M + V+... (3.4) 

with 



(k|V|k') 



1 



Wi + w 2 + w\ + w 



T (k\U\k') 



1 



1 



l^ 71 ") 4wwity 2 wi'u;2 



/r(k,k')<^ (3.5) 



In (|3.4|) the dots stand for higher order terms in a s and a and kinematic factors equal to 1 
on the energy shell have been neglected. 

From Eqs. (|3.3|) and ( p.7|) one obtains explicitly 



(k|t/|k'> 



\ 



(w 1 + w 2 )(w[ + w' 2 ) 

WiW2W l 1 w' 2 



Aa s 1 2 

3^2Q2^10?20 + q 



(Q ■ q) ; 

Q 2 



z 1 

+ 2Q2 k X k ' ' ( ai + G ^ + ^^[ g20 ^ 1 ' Q) ~ 1w{a 2 ' ^ + 

1 1/1 

+— <7i • (To H -<7i ■ CT 2 

6 4 V 3 



2Q 2 

(Q-q 1 )(Q-q 2 

Q 2 




(«i-Q)(« 2 -Q) 



+ 



(2vr) : 



rf 3 re iQ - r J inst (r,q,g 10 ,g 20 )} 



(3.6) 



with 



J mst (r,q,g 10 ,g 20 ) 



or 



9io + 920 



[qwvqw - q? + 9ioV to - q 2 :] + 



<?io<?20/ • Iq.t| . . IqTL 

•- -(arcsm h arcsm J J 



\ctr\ 



a 
r 



920 



?io - q| 



9io 920 

:(r x q-cTi + zg 10 (r -ai)) + 



(3.7) 



9io 



920 ~ qT 



r x q-(T 2 -«g 2 o(r-« 2 )) 



Here a*? denote the usual Dirac matrices 7° 7^, <r| the 4x4 Pauli matrices ( ^ 
obviously q= k 2 ^, Q = k-k' '" , + "''' 



9io 



l. ) and 
r ' i '-. Notice that, due to the terms in a*, such 



U is self adjoint only with reference to the undefined metric operator 7°7 2 . 

Due to (|3.5| ) the potential V can be obtained from U as given by (|3.6| - |3.7|) simply by the 

kinematic replacement 



(wi+W2){w' 1 +w' 2 ) 
wiW2w' 1 w' 2 



2\/ wiW2w' 1 w 2 



IV. DETERMINATION OF THE SPECTRUM 

In ref. || we have evaluated the eigenvalues of the operator Hcm for the potential V 
discussed above omitting the spin-orbit terms and including only the hyperfine splitting. 
The numerical procedure we have followed is very simple. It consists in solving first the 
eigenvalue equation for the static potential ( |1.1| ) by the Rayleigh-Ritz method |7| using 
the three-dimensional harmonic oscillator basis diagonalizing a 30 x 30 matrix. Then we 
have evaluated the quantities (V^I-HcmIVv) for the eigenfunctions tjj v obtained in the first 
step, choosing the scale parameter occurring in the basis in order to make minimum the 
ground state mass (V'is'I-^cmIV'is)- Notice that the determination of (Vvl^lVv) for the exact 
V is not trivial, since in general one should evaluate five-dimensional integrals of a highly 



6 



singular functions. For such reason we have used two different expansions for high and low 
transversal momentum (angular momentum), that allows to reduce to a three-dimensional 
integrals and treated the singularity with the method suggested in ||. 

The procedure we have followed for the determination of the eigenvalues of M 2 is es- 
sentially the same. Again we solve first the eigenvalue equation for Hqm with the static 
potential, and then we evaluate the quantities (ip u \M 2 \tfj u ). In this case the hyperfine split- 
ting is determined by the equation 

( 3 M nl ) 2 - CM nl ) 2 = 



97T JO JO V WiWo V W\Wn J -I 



which is more complicated than the corresponding equation in the case of if cm @ • 

Both the new results based on M 2 (crosses in figs.[l|-|3| and dashed lines in fig.|j) an d the 
old ones based on Hqm (circlets in figs.Q-(^ and dotted lines in fig.|]) are reported in figures 
for the parameters discussed in the introduction. For the / > cases masses represent the 
center of mass of the multiplets. In both cases the agreement with the data is on the whole 
good, not only for bottonium and charmonium (as in ordinary potential models), but also 
the light-light and light-heavy systems. Notice that the quadratic formulation seems to give 
a better low angular momentum light-light spectrum, while perhaps the linear formulation 
gives better Regge trajectories. The Regge trajectories in the quadratic case can be improved 
rising the value of o at the price, however, of making more difficult the fitting of the light- 
heavy states (perhaps the MAL model is too naive). As we mentioned the only serious 
disagreement remains that of the light pseudoscalar mesons. 
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FIG. 1. Heavy- heavy quarkonium spectra 
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FIG. 2. Light-light quarkonium spectra. Question marks indicate states whose assignment to a 
multiplet is not obvious. Position of rjs is derived from 77 and 7/, by standard mixture assumption 
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FIG. 3. light-heavy quarkonium spectra 
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